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We study the interplay between forgetful and memory-keeping evolution enforced on a two-level 
system by a multi-spin environment whose elements are coupled to local bosonic baths. Contrarily to 
the expectation that any non-Markovian effect would be buried by the forgetful mechanism induced 
by the spin-bath coupling, one can actually induce a full Markovian-to-non-Markovian transition 
of the two-level system's dynamics, controllable by parameters such as the mismatch between the 
energy of the two- level system and of the spin environment. For a symmetric coupling, the amount 
of non-Markovianity surprisingly grows with the number of decoherence channels. 
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I. INTRODUCTION 

The understanding of the implications of non- 
Markovianity and the reasons for its occurrence are still 
largely elusive. Yet, they are stimulating a growing in- 
terest in light of their potential impact on many disci- 
plines, from quantum information and nano-technology 
up to quantum biology [IH3]. An important contribu- 
tion to this quest came from the formulation of quanti- 
tative measures of the degree of non-Markovianity of a 
process [4]-[7]. In general, these tools address different 
features of non-Markovianity, from the lack of divisibility 
of a map [5] to the ability of the environment to recipro- 
cate the information transfer from the system. This pro- 
cess occurs unidirectionally in a Markovian dynamics [6] , 
while the re-focusing of information on the system is the 
signature of memory effects, as verified in all-optical set- 
ups [8]. The handiness of such instruments has recently 
triggered the analysis of non-Markovianity in quantum 
many-body systems such as quantum spin chains [9] or 
impurity-embedded ultra-cold atomic systems [10 and 
in excitation-transfer processes in photosynthetic com- 
plexes [3 . While these studies relate non-Markovian fea- 
tures to the critical behavior of a quantum many-body 
system [9j [11] , they also provide a promising arena where 
the roots for non-Markovianity can be researched in phys- 
ically motivated contexts. 

In this paper we explore the competition between two 
profoundly different mechanisms in a simple open quan- 
tum model that is relevant for the physics of nitrogen- 
vacancy centers in diamonds [12] and molecular nano- 
magnets [13 . Specifically, we address the interplay be- 
tween the dynamics induced on a two-level system by 
its coherent interaction with other (environmental) spins, 
and the Markovian process describing the relaxation of 
the latter. One would expect that, when such memory- 
less dissipative coupling determines the shortest dynami- 



cal timescale of the system, Markovianity should emerge 
preponderantly, especially as the number of environmen- 
tal spins increases. Indeed, one could imagine that a 
sort of "Markovianity-mixing" property would hold as a 
result of the increasing difficulty to re-build the coher- 
ence of the system when many decoherence channels are 
open. Quite strikingly, we show that this is not gener- 
ally true. In order to do this using a physically relevant 
model, general enough to encompass the unexpected fea- 
tures that we would like to highlight, we consider a spin- 
star configuration whose peripheral sites are coupled to 
rigid boson environments, assumed to induce a memory- 
less dissipative dynamics. While certainly not exhaust- 
ing the possible scenarios that can be tackled, our choice 
is illustrative since the degree of non-Markovianity (as 
defined in Ref. [6 ) can actually increase with the num- 
ber of peripheral spins, while stronger interactions with 
the boson baths only affect its rate of growth. The fea- 
tures of the system at hand are quite complex and a rich 
non-Markovianity phase diagram emerges, spanning de- 
grees of memory-keeping effects all the way down to zero 
values. This can be exploited to qualitatively modify 
the character of the dynamics by engineering its features 
via accessible control parameters such as the detuning 
between the central and the outer spins. In turn, this 
opens up the possibility to implement qubit-state prepa- 
ration protocols in an open-system scenario that exploits 
non-Markovinity, along the lines of Refs. [17 and be- 
yond the well-established Markovian dissipative frame- 
work [H[T9]. 

In the following, we first present the model and its 
solution in the simplest terms in Sec. [II] while the 
microscopic description and more sophisticated solution 
method are presented in the Appendices. We then pro- 
ceed to the analysis of the non-Markovianity of the dy- 
namics in Sec. |III|an d Sec. [TV) Some concluding remarks 
are given in Sec. |V| 
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II. THE MODEL AND ITS SOLUTION 

The physical set-up that we describe is sketched in Fig. 
[I]a, which shows a central spin (labelled 0) coupled to N 
outer spins, with bonds along the branch of a star. Each 
environmental spin is further coupled to a local boson 
reservoir. The evolution of the central spin is ruled by 
the master equation 

N 

dtpo® =Tr s {-i[H,p(t)] + J2 l >Mt)}} (1) 

i=i 

with p(t) the density matrix of the whole system. Each 
Lindblad superoperator Cj describes local dissipation at 
temperature T (the same for all the baths) as [14] 



= 7(« + " i&t&i, Pi/2) 

+ 7n(<r+p<7~ - {&Ta+, p }/2), 



(2) 



where 7 describes the effective coupling of each exter- 
nal spin to its thermal reservoir, populated by n = 
(eP^i — excitations (fi = l/k b T, where k b is the 
Boltzmann constant). In what follows, we will consider 
the peripheral spins to be initially prepared in <S>^ =1 | — ) ■. 



environment 
of spinj 




FIG. 1. (color on-line): (a) Scheme of the system: The central 
spin interacts with N peripheral spins, each affected by its 
own local environment, (b) Evolution of states \+) x (red tra- 
jectory) and \—) x (blue one) for a star with N = 4 peripheral 
sites. The Bloch spheres in the left (right) column corre- 
spond to the isotropic (anisotropic) spin-spin coupling. The 
top (bottom) row is for the resonant (off-resonant at A = 1/2) 
case with 7/7=1 and T = 0. In the isotropic cases, the final 
state of spin is pure, while for A = ±1 it is mixed. A/0 
prevents the intersections of the trajectories, which are the 
dynamical points at which the trace distance is strictly null. 

To solve the master equation, we use the damping ba- 
sis [16] made out of tensor products of eigenoperators of 
Cj. In this basis, the density matrix of the system reads 



where /}] = (£ 



<j?)/2, A? = *?/2, m = ^ and 



ftj — <r ■ are right eigenoperators of Cj with eigenval- 



ues Xj 



0,A 



-(2n + 1). The set of op- 



erators {O m } is composed of the tensor product of TV 
damping-basis elements, one for each peripheral spin. 
Due to the symmetry of the Hamiltonian, if O r and O s 
consist of the same elements of the damping bases (al- 
though differing for their order), the respective coeffi- 
cients must satisfy c nr = c ns . This simple observation al- 
lows us to reduce the number of relevant operators from 
4 N to N ?z N 3 . With the help of the single-spin dual 
damping basis {p™} (n = 1,..,4), made of left eigen- 
operators of Cj's, and using the orthogonality condition 
Tr[/^,/2f'] = S k yS jU we find 



N 



s(t) = c nm (t)M 

r 



(4) 



With M n mrs= iTl{(p^O s ) [H, &Q®dm} }+Am ^rn ^ms 

and A m = Y^=i x f- B ^ callin & = eM \ the state of 
the spin star at time t is 



p(t)= ^ JCnmrs(t)c r s(0)a^6 s . 



(5) 



r,s,n,m 



Tracing over the degrees of freedom of the peripheral 
spins, we find 



Po(t) = Yl ( X^mrlW C rl (0)a° 



C2l(t) 



V 



2n 2 

C3l(t) 



C3l(*) 

(1 + n) _ C2i(t) 
1 + 2n 2 



(6) 



This gives the exact solution for the dynamics of the cen- 
tral spin, valid for any N once the expressions for c rs (t) 
are taken. With this at hand, in the next section we eval- 
uate the amount non-Markovianity of the time evolution. 



III. NON-MARKOVIANITY 

To quantify the degree of non-Markovianity of the dy- 
namical evolution of the central spin described in Eq. (|6| , 
we employ the measure put forward in Ref. [6], which 
is based on the idea that memory effects can be char- 
acterized by the information flowing out of the open 
system and quantified in terms of the trace distance 
D[p OA (t),po^(t)] = Tr|p ,i(*) - Po, 2WI/2 between any 
two of its states Poj(t) (j = 1,2). The trace distance 
quantifies the distinguishability of two states and leads 
to measure non-Markovianity as 



p(t) = ^ C nm(*)£o ® On 



(3) 



n=l m=l 



Af — max 

Po,i(0) 



d t D[p ,i(t),p 0: 2(t)}, 



(7) 
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where Q + is the union of the intervals where d t D>0. To 
provide a general assessment of the dynamics of spin 0, 
we consider the coupling with the external spins to be 
described by the anisotropic XY model 



N 



H S = J^J(1 + \)a%a* + (1 - AKaJ] 

= 1 



(8) 



where A is an anisotropy parameter and J is the spin- 
spin coupling strength. For isotropic coupling (A=0) and 
zero temperature, we obtain a simple scaling law [15]: for 
any N>1 po(t) is obtained from the expression valid for 
N = 1 with the re-definition J — >> Jy/~N. This enables 
the analytic optimization over the input states entering 
M. By calling p^ l ~(k\po yi \l) , we have 



D[p 0!l (t),p 0:2 (t)}=V5p^(t)\g (t)\^+Sp^(t)\g (t)\ (9) 



where 5p kl (t) 
and we have 
exp[-i(G + iA)i][(G 



pUt) 



<iW - 

introduced gw(t) = 
iA) sinh(zt) + z cosh(zt)]/2z, 
s = ^(G + iA) 2 - J 2 N/2,G = 7 (n + 1/2) and 
the energy mismatch A = e — cq between the cen- 
tral and outer spins. The maximum in Eq. ^ is 
achieved for the pure states po,i = IV^o^il w ^h 
|^) = cos(^/2)|-) + e^sin(^/2)|+) . Here, (0*,&) 
are the angles that identify the respective Bloch vector. 
M is optimized by equatorial antipod al sta tes {i.e. states 
with #i 5 2 = 7r/2 and <p2 — tt — 0i). In VII, we provide an 



alternative analytic approach to the evolution of spin 
and the dependence of the trace distance on such angles. 

The trajectories described on the Bloch sphere by the 
evolved states are shown in Fig. [I] (b) [top row, left-most 
sphere] where we see that the states tend to intersect, giv- 
ing D = 0. For the states that optimize the mea- 
sure of non-Markovianity are those with (61,62) = (tt, 0) 
(the phases being immaterial) as shown in Fig.[2](b). In- 
terestingly, non-zero values of A hinder the intersections 
of the state trajectories [cf. Fig. [I] (b)]. However, this 
does not prevent the dynamics to become Markovian at 
proper working points, as we show later on. 

The evolution of spin can be characterized using J\f. 
When the peripheral spins are detached from their re- 
spective baths, any information seeded in the central site 
undergoes coherent oscillations from the center to the pe- 
riphery of the star and back. For A = 1 and peripheral 
spins prepared in 1/iV, the dynamics induced by RHsR 
with R = <Tq ®jLi cr V j is (strongly) non-Markovian at 
all times [6]. In our case, the interaction of the outer 
spins with their environments radically modifies this pic- 
ture. As an example, in Fig. [2] (a) we plot the trace 
distance for the optimal states at A = 0. We ramp 
up the spin-bath interaction strength 7, at set values 
of the intra-star coupling J, looking for the influences 
that an explicitly Markovian mechanism has on the de- 
gree of non-Markovianity that arises from the dynamical 
environment to which particle is exposed. We find a 
non-monotonic behavior of the trace distance that results 



in non-Markovianity. The quantitative features of D de- 
pend on the actual strength of the Markovian process: as 
7 increases, the revivals of the trace distance become less 
pronounced. As J\f depends on the number of temporal 
regions where d t D > 0, Fig. [2] (a) tells us that M de- 
creases as 7 increases, thus showing that, at resonance, a 
strong influence from the rigid environmental baths over 
the peripheral spins is sufficient to make the whole pro- 
cess Markovian. 

D(p ,i(t),po,2(t)) 
1.0 



... 7 /J = 0.5 

y/J= 1.0 
7 /J=1.5 




FIG. 2. (color on-line): (a) Trace distance between the opti- 
mal states po,j(t) (j = 1,2) for N = 6 peripheral spins with 
A = and 7/ J = 0.5 (dot-dashed line), 7/ J = 1 (dashed 
line) and 7/ J = 1.5 (solid line). As the relaxation time 
becomes shorter, the revivals of D(po,i (£), po,2(t)) are sup- 
pressed as a result of a reduction of information back-flow 
from the baths, (b) J\f against A for 7/ J = \/~N. The two 
lines correspond to #1,2 = tt/2 , (<pi, 02 ) = (0,7r) (solid blue 
curve) and (61,62) = (0, 7r) (dashed red curve), which are 
the optimal states in different detuning regions: J\f is the 
topmost curve in each region. There is a finite window of 
detunings (light-shadowed region marked as M) where J\f — 
[NM marks regions where J\f ^ 0]. Inset: J\f against A for 
7/ J = 0.5,0.75, 1, 1.25 and 1.5 (from top to bottom curve). 

This is expected as the excitations distributed to the 
peripheral spins by spin find the sink embodied by the 
baths. The reduced ability to feed back information sets 
J\T = 0. However, the general picture is more involved: 
it is sufficient to move to the off-resonant case to face a 
rather rich phase diagram of non-Markovianity. Fig. [2] 
(b) considers the case of coupling mechanisms such that 
7/ J = y/N and explores the effect that an energy mis- 
match between spin and the peripheral sites has on N '. 
We find two ranges of values of A for which J\f = 0, sym- 
metrically with respect to A = 0. In between and beyond 
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such regions, Af behaves quite distinctively: at resonance, 
the measure of non-Markovianity achieves a global maxi- 
mum (equatorial states realize the maximum upon which 
Af depends). For larger detunings, Af changes slowly with 
A (|±) being the optimal states). Clearly, the trend fol- 
lowed by Af also depends on 7/ J: small values of 7/ J 
push the dynamics towards strong non-Markovianity, re- 
gardless of A, as many coherent oscillations occur be- 
tween site and the periphery before the initial excita- 
tion is lost into the environments. At the same time, 
the range of detunings for which Af = increases with 7 
[cf. inset of Fig. [2] (b)]. However non-Markovianity per- 
sists, both on and off resonance, even when 7 becomes 
the largest parameter. This demonstrates an effective 
control of the degree of non-Markovianity of the dynam- 
ics undergone by spin 0, which can be tuned by both the 
energy mismatch between the outer spins and the central 
one, A, and the intra-star coupling strength., J. 




1 3 5 7 9 



FIG. 3. (color on-line): (a) Af against N and A for 7/ J = 
0.1 and n — 3. Differently from T = 0, except for a small 
range of values, the detuning has no effect on the character of 
the dynamics of spin 0. Strikingly, Af grows with N (almost 
linearly for N ^> 1). (b) Analytic behavior of Af versus N for 
A=0, A = 0, 7 = J at T — (n = 0). Inset: we present the 
case corresponding to A = — 1 [other parameters as in panel 
(b)]- 

Our discussions so far were restricted to the isotropic 
coupling at zero temperature, T = A = 0. When the pe- 
ripheral spins interact with baths populated by n thermal 
excitations, the Markovianity regions disappear. This is 
seen in Fig. [3] (a) where we show a typical case of the 
behavior of Af against A and N. The anisotropy of the 
intra-star coupling is crucial for the determination of the 
dynamics: for A^0 the pair of states that maximize Af 



changes with the number of peripheral spins. A numeri- 
cal search for the optimal states can be performed, lead- 
ing to quite surprising results concerning the scaling of 
Af with the size of the spin environment. Intuitively, 
one would conclude that, as N grows, the dynamics of 
spin will be pushed towards Markovianity. This is not 
the case: as shown in Fig. [3] (a), Af increases with N 
if A = 0, regardless of A. This shows that the non- 
Markovian character resists such Markovianity-enforcing 
mechanisms and, counter-intuitively, overcomes them. 
We have checked this behavior for the exact analytical 
expression obtained at A = A = [cf. Fig. [3] (b)]. The 
picture somehow changes for A 7^ 0: Af decreases with 
the growing dimension of the star. However, even for 
N ^> 1 the non-Markovian character is preserved and Af 
achieves a non-null quasi-asymptotic value. 



IV. TIME DEVELOPMENT OF 
NON-MARKOVIANITY 

The non-Markovianity measure gives an integral char- 
acterization of the dynamics. More details on the 
time dependence of the system-environment information- 
exchange process is obtained by considering the ratio 
of in-flowing to out-flowing information, up to a given 
value r of the evolution time. To this end, we de- 
fine K(t) = ^^y, where the in-flow [out-flow] A/" + (r) 
[Af~(r)] is defined as [minus] the integral of d t D, over 
the time intervals in which it is positive (negative), but 
only up to r. To evaluate these quantities explicitly, we 
chose as input states the same po,i that optimize the non- 
Markovianity measure Af = linv^oo A/" + (r). The ratio 
7Z(t) gives the fraction of the lost information that re- 
turns to the system within r, and its behavior is quite 
different in the various dynamical regimes that we have 
identified so far. In Fig. |4j 7Z(t) is shown for three val- 
ues of A corresponding to the three regions of Fig. [2] 
(b). The diverse evolutions of 1Z(r) signal qualitatively 
different dynamical behaviors of the system, depending 
on both the detuning and the anisotropy parameter. At 
short times, 7Z(r) is always zero (information has to flow 
out of the system before it can come back) , while its first 
peak is determined by the first revival of the trace dis- 
tance [see Fig. [2] (a)]. Then, its features become strongly 
dependent on A. At long times and at resonance, where 
a maximum of Af is found for A = 0, information oscil- 
lates between the star and spin Oand^(r) ^0 [cf. Fig. [I] 
(a)]. The overall dynamics is non-Markovian also for the 
case of Fig. [4] (c), where the time behavior of 7Z(r) is 
shown for a large detuning. In this case, however, 1Z{t) 
decays to zero at long times. Thus, the regions of non- 
Markovianity in Fig. [2] (b) correspond to different be- 
haviors: near resonance, a fraction of information comes 
back to the system, different input states remain dis- 
tinguishable even at long times and thus no equilibrium 
state is found. For large detunings, non-Markovianity 
is built up at short times, while different input states 
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FIG. 4. (color on-line): Ratio K(t) = J\f + (t) / J\f~ (t) versus 
A for a star of N = 8 sites at T = 0, with 7/ J = (left 
plots) and 7/ J = 1 (right plots) for three different values of 
the detunig: A/ J = for the plots (a) and (d), A/J = 0.7 
for (b) and (e), while A/J = 3 for (c) and (f). 



converge towards a long-time equilibrium. On the other 
hand, for intermediate values of the detuning [i.e. for 



A in the Markovianity region of Fig. [2] (b)] and A = 0, 
there is no back-flow. Even for A / 0, the fraction of 
information that comes back is quite small. The picture 
changes when J increases, the evolution becoming in- 
creasingly non-Markovian and the role of the anisotropy 
being fully reversed: A = implies a larger 7£(r), per- 
sisting for longer times at resonance. 

V. CONCLUDING REMARKS 

We have used a measure of non-Markovianity to show 
the possibility to control the dynamics of an open quan- 
tum system coupled to many independent decohering 
channels. We have highlighted the key role played by 
the detuning and the degree of anisotropy of the system- 
environment coupling: both can be used to explore a 
rich non-Markovianity phase diagram, where qualita- 
tively different scaling laws with the number of decoher- 
ence channels are found. The ability to switch from a 
Markovian to a non-Markovian regime by means of a lo- 
cal parameter could be used to prepare a quantum sys- 
tem in a desired state: indeed the Markovian character 
of processes can be employed for state engineering and 
information manipulation p~8j[19]. On other hand, while 
the formation of a steady entangled state is supported 
by non-Markovianity, a purely Markovian dynamics pro- 
duces separable steady states [T7] . 
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VI. APPENDIX A 

The total Hamiltonian of the spin-star system 
H=Ho+Hs+Hb consists of a few contributions. The 
first one is the system's free energy [we take units such 
that h=l throughout the paper] Hq = 52f=o 6 J '^j + 
52f=i 52k jbk,j-> describing the free evolution of 7V+1 
spin-1/2 particles [here, d s - is the s-Pauli matrix of spin j 
(s = x,y,z)], each with transition frequency ej between 
spin states and The second term in Hq de- 

scribes the energy of N sets of Mj harmonic modes (one 
set per peripheral spin of the star) with creation (anni- 



with 



hilation) operators b\ ■ (bkj) which satisfy the commu- 
tation relations [&fc,j, J = Skk'Sji- The central and 
peripheral spins are coupled by Hs, whose explicit form 
will be specified later on. Each peripheral spin interacts 
with its own bath as Hb = 52f=i Hbj, where 



with a- = (&j ± i&j)/2. We assume that the their lo- 
cal bath induce a Markovian dynamics of the peripheral 
spins and take uniform couplings, so that the evolution 
of the central spin is ruled by 

N 

dtPo(t) = Trs{-i[H,p(t)]+Y t £MW (11) 

3 = 1 



VII. APPENDIX B 

Here we provide an alternative solution to the dynam- 
ics of the system. In the interaction picture with respect 
to Hq the Schrodinger equation reads 

d t mt)) = -iH I (t)Mt)), (12) 

where the interaction Hamiltonian is given by 

N 

Hi{t) = ^E^w^ 7 W+^w<w) 



N 



J2T,(9k<(t)b k (t)+gtai(t)b{(t)) (13) 

j=i k 



af{t)=afe ±i ^ t (j = 0,...N) 
b k (t) = gkhe-™** 
b\{t) = gkble+^K 



(14) 



The operator N = 52j( a j + (52 k b\bk)j) counts the num- 
ber of excitations in the system and commutes with the 
total Hamiltonian H, so that any initial state of the form 



|*(0)) = (c |-) + c 1 (0)|+)°)|0> S |0) J 



N 



N 



^ Cj (o)|-)°|j) s |o) B + ^^ CfeJ (o)|-)°|o) s |k; 

j=l j=l k 

evolves after time t into the state 

l*W) = (co|->° + Cl (0)|+)°)|0) s |0) B - 



N 



N 



E^wi-) ij) s i ) B + EE c ^wi-) i°) s i k )f( 15 ) 

j=l j=l k 

where the state \0) S denotes the product state ®f =1 \— ) J 
and \}} S = (jj~\0) S for the sites on the star; \0} B is the 

vacuum state of all the reservoirs, and |k)^ = b£|0)j the 
state with one particle in mode k in the jth reservoir. 
The amplitude Co is constant in time because of 
^(t)|-)°|0) S |0) B =0. 



Substituting Eq. (15) into the Schrodinger equation (12) 
one finds 



N 



dt^J 



i(e —ej)t 



(Cj —UJkj)t 



i(ei—LJkj)t 



(16) 



We assume in the following that cj(0) — Ckj(0) — 0. This 
means that the two level systems on the star are in the 
|— ) state and that each environment is in the vacuum 
state initially. 

The total initial state is given by the product state 

l*(0)) = ( C o|->° + Cl (0)|+)°)|0) B |0) s . (17) 



Formally integrating Eq. ( |16| ) and substituting into 
Eq. ( [l6| ) one obtains the system for the amplitude 

ci(t),cj(t), 



N 



Cl(t) 



ijE c i(*) ei(e °" £i)t 

3=1 



f cj(ti)Ei^i ae<(£J " u * ,)(<_tl)(ft i 

JO 



(18) 
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We can define the kernels fj(t — t\) describing the two- 
point correlation function of each reservoir, which are the 
Fourier transform of the respective environmental spec- 
tral density 



f 3 (t-h) 



-U k j)(t — ti) 



(19) 



For the moment, we do not make any restrictive hypoth- 
esis on the form of fj , so that our results will be valid for 
an environment with a generic spectral density. In order 
to solve the system above it is convenient to pass in the 
Laplace domain: 



N 



sci[s] = c\ (0) — iJ Cj [s + z(ep 



(20) 



iJci [s — i(eo — cj 



}fj 



Co s \u s 



Solving the second of Eq. 20 respect to Cj [s], assuming 
that all the reservoirs are the same (fj(t) = f(t) V j), 
and substituting in the first we get 



ci \s 



:ci(0)- 



■ iA - f[s - i(e - e)] 



.is(e -e)-isf[s-i(e -e)] + J 2 N 



where A = (eo — e) (ej = e V j). 

To specify the model, but still retaining a general 
enough description, we consider a Lorentzian spectral 
density for each bath (which gives rise to an exponen- 
tially decaying correlation function): 



J(u>) 



1 



7A 2 



2tt (e j 



A 2 ' 



(21) 



Here 5 = Cj — uo c is the detuning of the center frequency 
of the bath uo c and the frequency of the two-level sys- 
tem ej , the parameter A defines the spectral width of the 
environment, which is associated with the reservoir cor- 
relation time by the relation tb = A -1 and the parameter 
7 is related to the relaxation time scale tr by the relation 

We will consider S = 0, and in this case we may dis- 
tinguish between the Markovian and the non-Markovian 
regimes (for the dynamics of the environmental spins 
themselves) using the ratio of 7 and A: 7 < | gives 
a Markovian regime and 7 > ^ corresponds to non- 
Markovian regime. 



Substituting in Eq. [21] and anti-transforming we have 
ci(t) = ci(0)G(t) with 



t*(t) = Ci(UJ 7 77 77 v , 

(ai - a 2 ){a 2 - a 3 ){a 1 - a 3 ) 

(22) 

where 5{ = a^— iA, j, k = 1, 2, 3 and for j<k. Here, a^s 
are the roots of the equation 



p(s) =s 3 + (2iA + X)s 2 



(23) 



(J 2 TV + A 2 + iAX + A 7 /2) 5 + J 2 N(iA + A) . 

Already at this point, it is evident how the only effects 
of increasing N is to redefine the coupling constant J. 

The solution of the Schrodinger equation of the total 
system with initial states of the form (17) lyes in the 



sector of the Hilbert space corresponding to zero or one 
excitations. 

We can construct the exact dynamical map describing 
the time-evolution of the reduced density matrix of the 
central spin which is given by 



p(t)=Tr s+ B{Mt))(*(t)\} 



Pu(t) pio(t) 
Poi(i) Poo(t) 



> (24) 



where pn(t) = (i\p(t)\j) for iJ = 0,1. Using Eq. ( |15| ) 
and Eq. (22) we find 



p 11 (t) = l-poo(t) = \c 1 (0)G(t)\ 2 , 
p W (t)=p* 01 (t) = c* oCl (0)G(t). 



(25) 
(26) 



The optimization of the initial states in Eq. ^ obtains 
the maximally possible non-Markovianity of a particular 
quantum evolution. 

In our case, the maximization is achieved by pure states, 
thus we choose as initial states for Eq. [TT] 

|*i(0)>= (cos(^)|-> + e 101 sin(0i)|+»|O> B |O) s (27) 

|* 2 (0)) = (cos(^ 2 )|-)+e^ sin((9 2 )|+»|0) B |0) s (28) 

With these states, the trace distance takes the form 

l -\G{t)\ V|G(t)| 2 (cos(^) - cos(^ 2 ))2 + ( S in(^) + sin(^ 2 ))2 

(29) 

where we used the fact that, since H is invariant under 
rotations along z-axis, the maximum is obtained for 

01 - 02 = TT- 



